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Abstract. 

Negative ion detachment in two-colour laser field is considered within the recent 

modification of Keldysh model which makes it quantitatively reliable. The general 

approach is illustrated by calculation of angular differential detachment rates, 

partial rates for particular ATD (Above Threshold Detachment) channels and total 

j\ detachment rates for H~ ion in bichromatic field with 1:2 frequency ratio. Both 

• »h . perturbative and strong field regimes are examined. Polar asymmetry and phase effects 

£>v are quantitatively characterized with some new features revealed. Phase effects are 

found to result in a huge anisotropy factor ~ 10 3 in the electron angular distribution 
in the perturbative regime. 



> 
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1. Introduction 

(N 

Photoionization of atoms in bichromatic laser field had received recently a considerable 

q^ \ attention both in theory (see, for instance, Baranova et al 1990, 1992, 1993, Baranova 

and Zel'dovich 1991, Szoke et al 1991, Anderson et al 1992, Schafer and Kulander 1992, 
Potvliege and Smith 1991, 1992a,b, 1994, Pazdzersky and Yurovsky 1994, Protopapas 
et al 1994, Veniard et al 1995, Pazdzersky and Usachenko 1997, Fifirig et al 1997) and 



experiment (Muller et al 1990, Ce Chen and Elliott 1990, Baranova et al 1991, 1992, 
Yin et al 1992). One of the principal points of interest seems to be dependence of the 
observables on the difference of field phases tp, i.e. the problem of phase control. Another 
important aspect is the angular (polar) asymmetry of photoionization rate. These effects 
are interrelated since polar asymmetry is ^-dependent and vanishes for some particular 
value of phase (see more detail below). The calculations have been carried out previously 
for ionization of hydrogen atom in two laser fields with a frequency ratio 1:2 (Schafer and 
Kulander 1992), 1:3 (Potvliege and Smith 1991) and 2:3 (Potvliege and Smith 1994). 
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Potvliege and Smith (1992) presented results for various frequency ratios and initial 
states. Different schemes were employed, but all of them implied numerically intensive 
work. 

For the multiphoton electron detachment from negative ions some analytical 
treatment exists (Baranova et al 1993, Pazdzersky and Yurovsky 1994, Pazdzersky 
and Usachenko 1997) being limited mostly to the case when one or both fields are 
weak. The presence of large number of parameters in the problem sometimes makes 
results of analytical studies not directly transparent. The case of fields with comparable 
(and large) intensities is also of interest bearing in mind both possible experimental 
realizations and the theoretical description of the transition between the multiphoton 
and tunneling regimes. 

The multiphoton electron detachment from negative ions presents unique situation 
when quantitatively reliable results can be obtained by analytical methods in a broad 
range of parameters characteristic to the problem. Indeed, it has been demonstrated 
recently (Gribakin and Kuchiev 1997a,b) that judicious modification of the Keldysh 
(1964) modelt ensures a very good quantitative agreement with results of numerically 
intensive developments, being much more simple. In many cases it also provides 
good agreement with available experimental data. It works unexpectedly well even 
for small number n of photons absorbed. In addition to numerous examples considered 
previously (Gribakin and Kuchiev 1997a,b), here we briefly comment on the most recent 
experiments by Zhao et al (1997) on non-resonant excess photon detachment of negative 
hydrogen ions. After absorption of two photons, the electron is ejected in superposition 
of S and D waves due to selection rules. The experiment demonstrates prevalence of D 
wave contribution (90% ± 10%). Our calculations give for D and S waves population 
86.2% and 13.6% respectively^ for experimental conditions (light frequency uo = 0.0428, 
field intensity / = 10 10 W/cm 2 . The elaborate numerical calculations by Telnov and 
Chu (1996) and by Nikolopoulos and Lambropoulos (1997) give for D wave population 
91% and 89% respectively. The difference between these values is almost the same as 
the difference between our result and that of Nikolopoulos and Lambropoulos (1997), 
all three theoretical predictions lying within experimental error bars. This, together 
with the cases considered earlier allows us to suggest that even for n = 2 an approach 
(Gribakin and Kuchiev 1997a,b) provides an accuracy comparable with that of the most 
elaborate numerical developments. 

The present paper extends approach of Gribakin and Kuchiev (1997a,b) to the 

f Subsequent development of this model was due to Perelomov et al (1966), Faisal (1973) and Rciss 
(1980); Perelomov and Popov (1967) were the first to consider multicolour process within this framework 
in terms of influence of higher harmonics on ionization probabilities. 

| The model shows also some population of G and higher partial waves. However, this unphysical 
effect proves to be less than 0.2% thus sustaining the model applicability. 



case of Dichromatic field. Its objective is to provide the scheme and some representative 
quantitative results for two-colour electron detachment from negative ions. In particular, 
the phase effects and the polar asymmetry are studied. The number of parameters in 
the problem is quite large and at the moment they cannot be fixed experimentally. 
Nevertheless it seems to be worthwhile to carry out some pivoting calculations in order 
to obtain insight into the possible magnitude of effects specific for negative ions in 
bichromatic fields. We consider angular differential detachment rates, heights of ATD 
(Above Threshold Detachment) peaks and total detachment rates. 

2. Scheme of calculations 

The previously developed scheme (Gribakin and Kuchiev 1997a,b) needs some 
modifications to incorporate bichromatic problem when electric field in the light wave 

p(t) = Fi coscjit + F 2 cos(uj 2 t + <p) (1) 

is a superposition of two harmonic components with the frequencies u>i, u> 2 and the 
amplitude vectors i*\, F 2 respectively; tp is the difference of field phases. Atomic units 
are used throughout the paper unless stated otherwise. 

We consider a case of commensurable field frequenciesf which implies that the 
common period T exists such that 

T=—M 1 = —M 2 (2) 

UJi UJ 2 

for some mutually simple integers Mi and M 2 . 

The exact expression for the differential transition rate dw\ is derived following 
Appendix A of the paper by Gribakin and Kuchiev (1997a) with the result 

dw x = 2tt J2 I A xe f | 2 KEx -E -e f )p x , (3) 

A ^f = ^J Q (Mt)\v(t)\Mt))dt. (4) 

Here i/jo(t) = exp(— iEot)(f>o describes an initial state for the time-independent 
Hamiltonian H , and ip\{t) is a quasienergy state for the total Hamiltonian H = 
H + V(t), which includes interaction with the periodic field V(t) = —er ■ F(t), 
V(t) = V(t + T): 

i^£ = [H + V(t)} Va, (5) 

Mt) = exp(-iE x t) 0a, M) = Ht + T ), (6) 

f The general treatment of incommensurable frequencies case was considered by Ho et al (1983), Delone 
et al (1984), Ho and Chu (1984), Manakov et al (1986), Potvliege and Smith (1992). 



E\ is the quasienergy, p\ is the density of states, f is the active electron vector. The 
energy ey absorbed by electromagnetic field could be presented as e/ = n\U\ + n<i u>2 
with some integers n\ and n,2- However this representation (i.e. the choice of n\ and 712) 
generally is non-unique which reflexes existence of different coherent interfering paths 
(with different number of absorbed and emitted photons of each frequency) leading to 
the same final state. 

If interaction of light wave with an electron is described in the dipole-length form, as 
presumed above, then a long range contribution to the matrix elements is emphasized. 
Therefore it is sufficient to employ an asymptotic form of the initial-state wave function 
(Gribakin and Kuchiev 1997a): 

o (f) w Ar"- 1 exp(-wr) Y lm (r) (r > 1), (7) 

where E = — |/t 2 , v = Z/k, Z is the charge of the atomic residual core (u — Z — 
for a negative ion), and r is the unit vector. The coefficients A are tabulated for many 
negative ions (Radzig and Smirnov 1985). 

The amplitude A\ e (4) is evaluated neglecting the influence of atomic potential 
on the photoelectron in the final state. Further on, the integral over time in (4) is 
calculated within the stationary phase approximation which presumes large magnitude 
of the classical action 

S(t) = \f (p+kt) 2 dt'-E,t, (8) 

where k t is the classical electron momentum due to the field 

k t = e f F(t') dt'. (9) 

The photoelectron translational momentum p plays a role of the quantum number A 
above; in particular, the quasienergy E\ = Ep, 

E *=¥ 2+ ik F > + i F > (10) 

includes contribution from the electron quiver energy due to the field. 

The result of calculations of the amplitude (4) could be written as a modification of 
the formula (25) in the paper by Gribakin and Kuchiev (1997a): 

4*, = -^^r(i + v/2) n> * Y.YUK) ^^ ■ ( n ) 

A corresponding expression for the detachment rate for the negative ion case [y = 0) 
reads: 

2 

dw ef 1 . |2 (27r) 2 



dfi p - " (27r) 2P l^ e/ l V 2 pA2 



E^(PV exp( ^ 



a/2^ 



(12) 



Here the subscript \x indicates that the function is calculated at the saddle point t = t^ 
which satisfies equation 

S%) = 0. (13) 

In the plane of the complex-valued time the saddle points t M lie symmetrically with 
respect to the real axis. Summation in (11) includes the points lying in the upper half 
plane (Imt M > 0) with < Ret M < T; p^ is a unit vector in the direction of p + k^. The 
magnitude of the final electron translational momentum p is governed by the energy 
conservation 

l -n 2 + E j; =e f) (14) 

which ensures that the momentum is real for open ATD channels. 
According to (8), (9), (1) we have 

S'(t) = 1 -(p + k t ) 2 + 1 -K 2 = 

1 e 2 e 2 

= o P 2 + ^ F ? sin2 w i* + TT^ F 2 sin 2 (w 2 t + <p) + 

6 6 

+p ■ Fi — sin uit + p ■ F 2 — sm(uj 2 t + ip) + 

- - e 2 1 

+Fi ■ F 2 sinc^it sin(a;2t + (/?) + - k 2 . (15) 

UJ\UJ 2 I 

The frequencies iO\ and uj 2 are integer multiples of the basic frequency u = 27r/T 

ui = Miu, ijj 2 = M 2 uj. (16) 

Assuming for definiteness that M 2 > Mi and introducing ( = exp(iut) we notice that 
the function 

V(C) = ( 2M *S'(() (17) 

is a polynomial of the power 4M 2 in (. This observation is of practical importance since 
equation (13) defining the saddle point is equivalent to 

no = o. (is) 

The efficient numerical procedures for finding roots of polynomials are available, and, 
in particular, one can be confident that all the roots are found. 

The practical calculations are conveniently carried out using the Mathematica 
program (Wolfram 1991). Starting from the expression for S'(t) one derives S(t) and 
S"(t) by analytical integration and differentiation respectively. The saddle points are 
found using Eq.(18), and the roots t M lying in the upper half plane are selected. Finite 
summation over fi in (11) or (12) completes the calculation. 
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The roots t^ and hence the photoionization amplitude (11) and rate (12) depend on 
the orientation of electron translational momentum p with respect to the field amplitudes 
F\ and F 2 . It is not difficult to consider the fields of various relative orientation and 
polarization, but for simplicity we limit our further calculations to linear polarized fields 
with F\ || F2. Then the differential photoionization rate depends only on the single angle 
9 between the vectors p and F\ || F2. 

3. Results 

Our calculations for H detachment are carried out for the parameters k = 0.2354, 
A = 0.75 (Radzig and Smirnov 1985). The frequencies ratio u>i/u 2 = 1 : 2 is considered. 
In this case the field (1) have zero mean value, but possesses polar asymmetry (i.e. 
asymmetry under inversion of the z axis directed along Fi \\ F 2 ) which could be 
conveniently characterized by the time-average value (Baranova et al 1990) 

(F 3 ) = -F?F 2 cos tp. (19) 

Presuming that Fi, F 2 > 0, from this expression one infers, for instance, that for 
the phase tp E [0, \n] the electric field effectively attains larger values in the positive- 
z direction than in the negative- ,2 one. This is illustrated, for example, by Fig.l in 
the paper by Schafer and Kulander (1992), or by Fig. 2 in the paper by Baranova et 
al (1993). Note that our definition of the phase tp is the same as in the papers by 
Baranova et al (1990), Muller et al (1990), Pazdzersky and Yurovsky (1995), but differs 
from that chosen by Schafer and Kulander (1992) who describe the electric filed as 
F(t) = Fi sinc^it + F 2 sin(w 2 t + tp KS ). Namely, the phases are related as <^ks — V 9 ~~ \ n - 

Although the formula (19) shows that the field has polar symmetry for tp — ^ir and 
the maximal polar asymmetry for tp = 0, quite paradoxically, the differential detachment 
rate (12) possesses polar symmetry for tp — (i.e. is invariant under transformation 
6 =^> 71 — 8), and is asymmetrical for other values of phase (see discussion by Baranova 
et al (1990), Schafer and Kulander (1992), Pazdzersky and Yurovsky (1995)). 

The other features of the phase effects are as follows. 

(i) The system Hamiltonian is a 27r-periodic function of the parameter tp. 

(ii) The system Hamiltonian is not changed by simultaneous transformation tp =^> — tp, 
9 =^> 71 — 0. The same applies to the differential detachment rate (12). 

(iii) The transformation tp =^> tx — tp leaves the Hamiltonian invariant only if t is replaced 
by-t 

As stressed by Baranova et al (1990), Baranova and Zeldovich (1991), Anderson 
et al (1992), Baranova et al (1993), the problem is invariant under the time inversion 



operation provided the final-state electron interaction with the atomic core is neglected. 
This is the case in the present model. Therefore our differential ionization rates are the 
same for tp and (ir — tp); hence it is sufficient for us to consider phases only from the 
interval tp G [0, |tt]. The calculations by Baranova et al (1990) within the perturbation 
theory and by Schafer and Kulander (1992) within the wave propagation technique 
took into account the final state electron-core interaction. Therefore they have found 
some deviations from the symmetry under tp =>- (ir — tp) transformation. However, for 
the negative ion photodetachment this effect could be anticipated to have only minor 
influence. 

At first we consider two fields with the frequencies u = 0.0043 and 2u and equal 
intensities I± — I 2 — 10 9 W/cm 2 . It is well known that the regime of detachment process 
is governed by the Keldysh parameter 7 = un/F (7 3> 1 for multiphoton detachment 
in perturbative regime; 7 <^C 1 for strong field, or tunneling regime). In the present case 
for the first field we have 71 = uj\KJF\ = 6, and for the second field 72 = 271, which 
corresponds to multiphoton regime. 

Fig.l (as well as Figs. 2-4 below) shows differential detachment rate as a function of 
the angle 0. The abscissas of the plots give the magnitude of 

1 1 dw ef 1 1 1 2 

2 2^d^0 = j2^ P \ Ap ^ ' (20) 

where the right hand side was calculated using the right hand side of the formula (12). 
The left hand side of Eq.(20) has the factor |. It means that the true detachment rate 
in case of H~ ion is twice larger than that given by Eq.(12). By this we account for 
the two possible spin states of residual H atom (i.e. for the presence of two equivalent 
electrons in H~). 

In Fig.l we show the differential detachment rate for the first and second ATD peaks 
which correspond to absorption of n = 7 and n = 8 photons of frequency u respectively. 
In Fig. 2 the same results are shown but for doubled value of the field amplitude F 2 . In 
Fig. 3 the amplitude F 2 is the same as in Fig.l, but the amplitude F\ is doubled. 

In all cases the angular distributions exhibit strong dependence on the field phase 
difference tp. This is well expected since the angular patterns are formed by interference 
of contributions coming from different complex- valued moments of time t M . For tp = 
the distribution is rather flat, with tp increasing it becomes more oscillatory. An 
interesting, and not obvious feature is that for tp — |7r the rate turns zero at the 
values of angle where it has minimum. 

In Fig. 4 we present the results for the same frequencies as before and equal field 
intensities I\ — I 2 — 10 11 W/cm 2 . Here the Keldysh parameter for the first field is 
71 = uj\KJF\ = 0.6. Bearing in mind the presence of the second field one can suppose 
that the situation corresponds to the onset of strong-field domain. The first open ATD 
channel corresponds to absorption of n — 18 photons of frequency 00. The patterns in 



differential rates become more oscillatory than in the weak filed case. 

The partial detachment rate for each ATD channel is obtained by integration of (12) 
over angles 

/dw e , f 9=0 dw e , , s 

-^T<%= / -^i-dcosO. (21) 

dilp J e=7T dcosV 

We present separately the result wi/ of integration over the upper half-space of 
the electron ejection (0 < 9 < |7r) and its counterpart w ef for the lower half-space 
(!"7r < 9 < 7r) . These magnitudes give a bulk characterization for the partial rate polar 
asymmetry. As discussed above, the polar asymmetry disappears (i.e. w^ = w N ) for 
p = 0. In the perturbative regime (Fig. 5) for the same conditions as in Fig.l we see 
that the bulk polar asymmetry parameter w^ /w N exceeds 10 3 even for the lowest ATD 
channel (N = 1) provided the phase (p is not too small (the open ATD channels are 
labeled by the number N = 1,2,... in order of increasing emitted electron energy). For 
higher ATD peaks the bulk asymmetry is even larger. 

The partial detachment rates integrated over all ejection angles w n = w^ + w N is 
shown in Fig.6 for three representative values of (p. Even for N = 1 the variation of the 
phase <p leads to the substantial change in the detachment rate described by a factor 4. 

In the tunneling regime the bulk polar asymmetry (Fig. 7) is not as prominent as in 
the perturbative regime. Nevertheless it is quite substantial. Note that both in Fig.5 
and Fig. 7 the electron emission in the upper half-space (0 < 9 < |7r) is more probable 
for all N, i.e. w^ > u>n m the interval of phases considered (<p G [0, \n\ ). As discussed 
above (see the property (ii)), the situation is reversed for <p e [ — ^-tt, 0] . For small value 
of the phase ip = |7r in Fig. 7 there is a clear tendency to swap the relation between 
wjy and w N for higher values of iV which is prevented by a kind of 'pseudocrossing'. 
For larger ip the partial rates w N are strongly suppressed when N increases. The phase 
effects in the partial rates wn are less significant in the tunneling regime (Fig.8). 

The total detachment rates are obtained by summation over all open ATD channels: 

w = ^2w N . (22) 

N 

The results for the total rates as well as w^ l ' u ^ = ^ N w^ u are presented in table 1. In 
the perturbative regime we see again strong bulk asymmetry (three orders of magnitude 
and more) if the phase difference p is not close to zero, and substantial variation of w 
with (p. Actually this result reiterates that for the partial rate with N — 1 since the 
latter gives a dominant contribution to the total rate in the perturbative regime. 

In the strong field regime the bulk polar asymmetry w^/w^ remains well 
manifested in the rate summed over all ATD channels. However, the total rate w is 
practically insensitive to the phase variation. The partial rates wn in Fig.8 exhibit 
some oscillatory structure as functions of N with position of extrema depending on p>. 
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This if- dependence almost completely disappears after summation over N as table 1 
shows. 

4. Conclusion 

As a summary, the approach of Gribakin and Kuchiev (1997a,b) provides convenient 
tool for investigating two-colour photodetachment of negative ions. The bichromatic 
electron detachment for H~ ion in the fields with 1:2 frequency ratio is examined in the 
perturbative and tunneling regimes. The polar asymmetry is found to be tremendously 
strong (~ 10 3 ) in the perturbative regime. Note that the asymmetry remains strong and 
keeps the same sign for all ATD for a wide range of phases < tp < it. This property 
makes this effect be convenient for experimental observation because it manifests itself 
very strongly even for a relatively poor resolution of the phase Stp ~ n/2 . It should 
be noted that via the recoil mechanism the predicted effect leads also to acceleration of 
the core thus creating anisotropic flux of neutral H atoms. 
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Table 1. Total rates w (summed over all ATD channels) for detachment of H~ ion in 
the bichromatic field with the frequencies u> = 0.0043 and 2w, equal intensities I\ = I2 
and some representative values of the phase difference ip. The detachment rates w^ u ' 
and w"' for electron ejection into the upper and lower half-spaces respectively are also 
shown. 



<p 


w 


w (u) 


W (D 


w 


w {u) 


W (D 




in units 10~ 21 a.u. 


in units 10 -6 a.u. 




I 1 = I 2 = 10 9 W/cm 2 


I 1 = I 2 = 10 n W/cm 2 





1.55 


0.773 


0.773 


162.1 


81.1 


81.1 


h 


2.12 


1.88 


0.240 


164.6 


99.8 


64.8 


¥ 


3.56 


3.50 


0.00566 


164.8 


113.0 


51.8 


¥ 


5.07 


5.06 


0.00112 


165.1 


123.8 


41.4 


¥ 


5.71 


5.71 


0.00040 


164.9 


128.0 


37.0 
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Figure captions 



Figure 1. Detachment of H~ ion in bichromatic field with the frequencies ui = 0.0043 
and 2ui and equal intensities I\ — I2 — 10 9 W/cm 2 . Differential detachment rate (see 
Eq.(20)) (in units 10 _12 a.u.) as a function of the angle 9 is shown for the first (a, 
absorption of n = 7 photons of frequency 10) and second (b, n = 8) ATD peaks and 
various values of the field phase difference: solid curve - cp = 0; short-dashed curve 
- f = tt7t; dot-dashed curve - ip = in; dotted curve - if — |7r; long-dashed curve - 
f= |?r. 



Figure 2. Same as in Fig.l, but for unequal field intensities I\ = 10 9 W/cm 2 , 
I\ = 4 x 10 9 W/cm 2 . The differential detachment rates is shown in units 10~ 10 a.u. 



Figure 3. Same as in Fig.l, but for unequal field intensities I\ = 4 x 10 9 W/cm 2 , 
I\ = 10 9 W/cm 2 . The differential detachment rate is shown in units 10 _10 a.u. 



Figure 4. Same as in Fig.l, but in the strong field regime: I\ = I2 = 10 n W /cm 2 . 
The detachment rate for the first (a, absorption of n = 18 photons of frequency lo) 
and second (b, n = 19) ATD peaks is shown in units 10~ 6 a.u. 
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Figure 5. Partial detachment rates for various ATD channels for H ion in 
Dichromatic field with the same parameters as in Fig.l (perturbative regime) and 
various values of the field phase difference <p. N labels ATD peaks with the lowest 
N = 1 peak corresponding to absorption of 7 photons of frequency u> = 0.0043. The 
rate w N of electron emission in the upper half-space is shown by circles, its counterpart 



w 



for the lower half-space is depicted by triangles. The plot for ip = g7r additionally 
includes the rate for ip = (crosses) when emission is polar symmetrical (w^' = w^J ) . 
The symbols are joined by lines to help the eye. 



Figure 6. Same as in Fig. 5, but for the detachment rate integrated over all angles 
= Wff + w N . The results for three v; 
0; blocks — <p = \tt, triangles - <p = |tt- 



wn — wjy + w N ' . The results for three values of the phase p are shown: circles 



Figure 7. Same as in Fig. 5, but for the field parameters chosen as in Fig. 4 (strong 
field regime). The lowest N = 1 peak corresponds to absorption of 18 photons of 
frequency u = 0.0043. The detachment rate is shown in units 10 -6 a.u. 



Figure 8. Same as in Fig. 6, but for the field parameters chosen as in Fig. 7 (strong 
field regime). The detachment rate is shown in units 10 -6 a.u. 
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